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Abstract: We prove nonexistence of a nontrivial integral that is polynomial in momenta 
of degree <6 for the Zipoy-Voorhees spacetime with the parameter 5 = 2. 

1. Introduction and the result 



We consider the following 4-dimensional metric on an open set U C 



p4. 



'1.1^ 



+ (x2-l)(l-J/2)rf^2 



X + 1 



dV 



This metric was introduced in [3l H]. It is Ricci-flat, has Lorenz signature, and two com- 
muting Killing vector fields: dt {ds^ is static) and ids^ is axially-symmetric) . 

Our interest in studying this metric is motivated by the resent papers of J. Brink [5l |6l [71 E] • 
She constructed numerically the geodesies of this metric and found out that they behave 
quite unusual. As it is clearly seen on figure 6 in [6], the behavior of the geodesies is in no 
way chaotic and is typical for the behavior of trajectories of a Liouville-integrable system. 
Later, in a private conversation. Brink demonstrated us that this 'integrable' behavior 
happens for all geodesies. 

Recall that a function I : T*U — M is called an integral for the geodesic flow of the metric g 
if it Poisson-commutes with the kinetic energy H = \g^^PiPj- Denoting by { , } the canoni- 
cal Poisson bracket this condition writes as {/, H} = 0, and it means that the function / is 
constant on the trajectories of the Hamiltonian system. In other words, / is the conserved 
quantity of the Hamiltonian H. A metric g is Liouville-integrable (for dim U = 4) if, in 
addition to the Hamiltonian, there exist 3 Poisson-commuting functionally independent 
integrals for the geodesic flow. Recall that integrals are functionally independent if their 
differentials are linearly independent almost everywhere. 

Metric (11. ip already has two additional commutative integrals, namely the momenta p^ 
and pt] the numerical behavior of the geodesies is an indication of the existence of a fourth 
integral, commuting with pt and p^. 

In the papers [71 [8], Brink suggested to look for an integral of the geodesic flow that 
is a homogeneous polynomial in momenta (with coefficients depending on the position). 
There are two reasons for this choice of an ansatz for the integral. The first one is that 
for all physically interesting systems the known integrals are such. The second reason is 
more mathematical: since the metric is real-analytic, it is expected that the integral is 
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real-analytic as well. Now, as it was known already to Whittaker [TT], the existence of an 
integral analytic in momenta implies the existence of an integral that is a homogeneous 
polynomial in momenta (More precisely, ii F = X^^g is an integral, where Fk is a 
homogeneous in momenta of degree k, then every Fk is an integral. If F commutes with 
Pt and p^, then every Fk commutes with pt and p^). Thus we can restrict to the integrals, 
which are polynomials of degree k in momenta. They are essentially the same objects as 
Killing (0, A;)-tensors. 

The main statement of our paper is that, contrary to numerical observations, no nontrivial 
integral of the expected low degree exists. The 'trivial' integrals for the geodesic flow of 
the metric g are H, Ii = prp, I2 = Pt and any (polynomial) function of them. 

Theorem. There exists no smooth function I : T*U — which is a polynomial in 
momenta of degree < 6, such that H , Ii, I2, I are functionally independent and Poisson- 
commute. 

This is a rigorous mathematical statement and its proof is also rigorous, though it is heavily 
based on computer algebra calculations. In the next section we explain the mathematical 
foundations and the details of these computations. The calculations are too complicated 
to be presented here in all the details, but the idea behind them is rather simple and can 
be easily realized in most computer algebra packages (the Maple worksheet is available 
from the authors). 

These computations can be extended to decide existence of higher degree integrals - the 
algorithm is unchanged (in fact, we have also verified nonexistence of the integral of degree 
7, but the hardware capacity limited us in degree 8). We will give some further comments 
in the conclusion. 

2. Nonexistence of the integrals 

In this section we demonstrate nonexistence of a nontrivial integral of degree 6. This in 
turn implies nonexistence of a nontrivial integral of smaller degree. Indeed if such an 
integral / existed in degree /c < 6, then Pi~^I would be a nontrivial integral of degree 6. 

2.1. The general idea behind the calculations. The condition that the function 

(2.1) /= I^3kmPlP^yP>T 

commutes with H is equivalent to a linear system S of PDE of the first order on the 
coefficients lijkm'-, the latter are smooth functions on U. Using the assumption that / 
commutes with p^ and pt, we see that the unknowns lijkm are functions of x and y only. 
More details on the differential equations in S will be given in the next subsection. 

Let us now take this system S and differentiate its equations with respect to the variables 
X and y, denoting the results by Sx and Sy respectively. The system S^^'^ = S U U Sy 
is called the first prolongation of 5". Every smooth solution of S is of course a solution of 
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S^^\ Next we consider the second prolongation S^"^^ = S^^^ U S^i U 5*2,.^ U Sy2 and so on. 
The n-th prolongation is S'^") = IJfc+m<n S.j.kym. 

The system 5" has finite type if for some £ one can solve the equations S''-^^ with respect 
to all highest derivatives of the unknown functions (in our case with respect to all jets 
of order £ + 1). It is known [TUl IH] that our system is of finite type. To see this let us 
notice that in the geodesic coordinates the metric is fiat to the second order at the given 
point. Since the coefficients of the top derivatives of the unknowns in the system S (and 
its prolongation S**^"-*) depend only on the zeroth and first derivatives of if, its symbolic 
behavior at any point is the same as for the fiat metric. In other words, the symbols of 
this system are isomorphic at all points and achieve finite type at the same level t. 

The calculation for the fiat metric is not difficult, and we obtain i = k ioi the system 
describing integrals of degree k. In particular, for our S we need to prolong £ = 6 times to 
achieve finite type (a calculation supporting this claim will be shown in §2.2p . 

For every linear system of finite type, the space of solutions is a finite dimensional vector 
space. Indeed, by the classical argument, the solutions are given as integral surfaces of the 
Cartan distribution, which in our case has rank 2 (on the equation-manifold given in the 
space of {i + l)-jets by S^^^). Thus a solution / of 5 is uniquely determined by the values 
of its derivatives up to order £ + 1 (here 1 is the order of S) at any fixed point (xo,yo)- 
This follows from the observation that restriction of / to a curve 7 C M^(a;, y) through the 
point {xq, yo) reduces 5" to a system of ODE in the Euler form. 

Take n > i and consider the prolonged system S^'^\ Treating the derivatives of lijkm at 
{x, y) up to the order n + 1 as independent variables u (called jets) allows us to write this 
linear system in the form Au = 0. Because coefficients of the Zipoy-Voorhees metric ds'^ 
are algebraic (rational) functions, the entries of the matrix A are also algebraic as functions 
of the point {x,y). 

If at a certain point {xq, yo) the rank of the matrix A equals to the dimension of u, the only 
solution is trivial u = 0. From the finite type condition it follows then that the solution u 
is identically zero. 

As we explain in the next section, our PDE system 5" decouples into two linear subsystems 
which have to be solved independently. For one of them Sodd the situation is as described 
above: for the 5*^ prolongation the linear system Au = has only the trivial solution. For 
the second subsystem Sgyen the situation is more complicated, since it has 16-dimensional 
space of solutions. These solutions correspond to the following 'trivial' integrals: 



2 4 6 

(2.2) = aH^ + H'-Y, mil'" + ^^-^1^2"' + E ^2 

i=0 j=0 k=0 



where a, Pi, Vj^rjk are arbitrary constants. 
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Thus for Seven the matrix A has a 16-dimensional kernel, whence rk{A) < dim(ii) — 16. If 
rk{A) = dim[u) — 16, then every solution u of the equations Au = is 'trivial' because the 
system is of finite type and n > i. In other words, u has form fl2.2p globally. 

How do we check that the finite type level £ does not exceed the number n we consider? 
If for a certain n we have rk{A) = dim(ti), the condition ^ < n is fulfilled automatically 
since we can express from the equations in S^""^ all derivatives of the unknown functions. 
This finishes the case Sodd- In fact in this case i = 5 (this is a simpler system than Seven 
and it achieves finite type earlier) and we go to tt, = 5 in our calculations. 

For Seven for n = Q we will have dim(u) — rk{A) = 16, which equals to the dimension of 
the space of trivial integrals of degree 6. Since £ = 6 by the general theory, we conclude 
that every integral must belong to the family (12. 2p . which finishes our argument. 

The size of the matrix A in our calculations is quite big, and it is not possible to handle 
it by hand. We use the symbolic software Maple. This is possible as the idea described 
above can be realized algorithmically. 

Actually we perform prolongations (symbolic differentiations) and operate with polynomial 
or rational functions. The coefficients of the latter are rational numbers and the point we 
substitute is also rational - we choose {xq^Uq) = (|,2). Then we calculate ranks of the 
matrices with rational numerical entries, which is done via the Gauss method. All of these 
calculations are exact (no approximations), and so the result has to be considered as a 
computer assisted mathematically rigorous proof. 



2.2. Details of calculations. We look for an integral I of form (12. ip with coefficients 
lijkm being smooth functions on U. Since = {I,h} = 0, these coefficients do not 

depend on 0, t and are smooth functions on the open domain in M^(x,?/) (given by the 
conditions x, y, ^ ±1, x ^ ±j/). 

Since the coefficients of H also do not depend on 0, t, the commutation of H and / writes 
as 

(2 3) [H I] - I ^^^^ -Q 

dx dpx dx dpx dy dpy dy dpy 

Substituting formulae for H and / into this expression we get a polynomial in momenta of 
degree 7. Its vanishing is equivalent to vanishing of all its coefficients, and there are 120 of 
them. These coefficients are rational functions of (x, y) with integer coefficients and they 
are linear in lijkm and their first derivatives. The number of unknowns lijkm = hjkm{x,y) 
is 84, and thus (12. 3 p gives a first order linear overdetermined system S of PDE. 

Notice that the Hamiltonian H is even in the total degree by the variables Pcf,,pt (i.e. 
it does not contain quadratic terms PxPt and similar). Then writing / = lodd + hven, 
where the index refers to the parity of A; + m in (12. ip . we observe the splitting {H, 1} = 
{H, lodd} + {H, leven} and the summands in the last expression preserve the parity (e.g. 
{H, lodd} has odd total degree by the variables p,p,Pt)- 
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Thus our system of equations decouples: S = Sodd^ Seven, Sodd^ Seven = 0- The subsystem 
Sodd is a system of hnear PDE on lijkm with k + m odd, and the subsystem Seven is a 
system of hnear PDE on hjkm with k + m even. Both Sodd and Seven are overdetermined. 

In fact Sodd consists of 60 equations on 40 unknown functions, while Seven consists of 60 
equations on 44 unknowns. Our goal for the subsystem 5*0^^ is to show that its certain 
prolongation at some point {xq, yo) satisfies dim{u)—rk{A) = 0. Our goal for the subsystem 
Seven is to show that its certain prolongation at (xq, yo) satisfies dim(n) — rk{A) = 16. We 
follow the scheme described in §2.11 and let {xq, ?/o) = (|, 2). 

Now we put the results of our calculations into the table. The first row in the table is 
the number of equations in S"^^], the second is the number of the unknowns lijkm (with 
k + m odd) and their derivatives by x, y up to order n + 1. The third row is the rank of 
the corresponding matrix A of the size of eqn) x (dim(ti)). 

The number of nontrivial integrals is the minimum of the quantity A = dim(u) — rk(A). 



n 





1 


2 


3 


4 


5 


# of eqn 


60 


180 


360 


600 


900 


1680 


dim(n) 


120 


240 


400 


600 


840 


1440 


rk{A) 


60 


180 


360 


590 


838 


1440 


A 


60 


60 


40 


10 


2 






We see that the 5* prolongation is enough to prove that the system Sodd has only trivial 
solutions. 

Next comes the table for Seven- The meaning of the rows is similar. The number of 
nontrivial integrals is the minimum of the quantity A = dim(u) — rk(74) — 16. 



n 





1 


2 


3 


4 


5 


6 


# of eqn 


60 


180 


360 


600 


900 


1260 


1680 


dim(n) 


132 


264 


440 


660 


924 


1232 


1584 


ik{A) 


60 


180 


360 


600 


888 


1215 


1568 


A 


56 


68 


64 


44 


20 


1 






We see that the 6* prolongation is enough to prove that the system Seven has only 'trivial' 
solutions, i.e. every integrals leven is of the form (12.21) . 

Finally we demonstrate how to calculate the number i symbolically. As explained in §2.11 
it is enough to compute it for the fiat metric. We form equation S similar to the above, but 
now we count only the equations of order n + 1 in S^"'\ and we separate top-order terms. 
The latter are the derivatives of lijkm of order n + 1 denoted by v, and the subsystem 
writes as a linear inhomogeneous equation Bv = w, where w combines derivatives of order 
< n. When A = dim(f ) — rk{B) = the corresponding n = i. Here is the table for Seven- 
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n 





1 


2 


3 


4 


5 


6 


^ of eqn 


60 


113 


166 


219 


272 


325 


378 


dim(f ) 


88 


132 


176 


220 


264 


308 


352 


rk(5) 


60 


113 


166 


214 


262 


307 


352 


A 


28 


19 


10 


6 


2 


1 






Thus £ = 6. Similar calculations for Sodd yield £ = 5. This finishes the proof of the theorem. 

3. Conclusion 

The question of existence of an additional integral is crucial for understanding of static 
axially-symmetric Ricci-flat metrics. The most famous spacetimes in the general relativity, 
given by the Schwarzschild and Kerr metrics, can be effectively studied because their 
geodesic flows admit an additional quadratic integral. This integral allows to describe 
and to control the behavior of the geodesies, and it also helps solving the wave and other 
physically-relevant equations. 

The family of Zipoy-Voorhees metrics with a real parameter S is given by the formula 

Transformation (x, 6) ^ {—x, —6) is the symmetry of this family, so we let 6 > 0. Pa- 
rameter 6 = corresponds to the flat metric, while for 6 = 1 we get the Schwarzschild 
metric. These are the only two cases from the family that admit Killing vector fields. 
For the Schwarzschild spacetime the symmetry algebra is in fact non-commutative, but its 
nontrivial Casimir function is an additional integral that yields Liouville integrability. 

We have checked that no other metric from the family admits a quadratic integral. The 
search for higher degree integrals we restricted to the next interesting parameter 6 = 2, 
corresponding to the metric (11.11) . We hoped that the integral of degree 4 exists; some 
other groups of mathematicians have also tried to find the missing integral by looking for 
a lucky ansatz. But, unfortunately, we have not found it and moreover we have proved it 
cannot exist in degree < 7. Degree 7 we have investigated with other methods (not the 
one presented in the paper), which also yield nonexistence. 

We stopped at degree 8 by purely technical reasons - even though the operational memory 
was not suffering, the computation time got too long. Actually, Maple obtains the matrix 
A quite fast; the hard part of the calculations is to compute the rank of the matrix A, 
which for the integrals of degree 6, has size 1680 x 1584. Our standard PC required few 
days for this. But we hope that a faster computer and more specialized symbolic software 
(which handles rational numbers in exact manner) could advance further. 

Let us remark that the Zipoy-Voorhees metrics can be included in the bigger family of 
Manko-Novikov metrics, which are also given by explicit formulas. Every axially-symmetric 
stationary Ricci-flat metric can be written with the help of the Ernst equations. Having 



NONEXISTENCE OF THE 6TH DEGREE INTEGRAL FOR ZIPOY-VOORHEES METRIC 7 

explicit formula for the integral for the Zipoy-Voorhees spacetime would suggest a per- 
turbaton of this formula to integrate other stationary axially-symmetric Ricci-flat metrics. 
J. Brinck calculated numerically also the geodesies of certain Manko-Novikov metrics and 
indicated numerically an integrable behavior. 

Existence of an integral of high degree is not completely impossible, as recently some non- 
trivial superintegrable systems were found. Its search can follow the proposed algorithm. 
One can further simplify it as it is apparent that the matrix A contains a lot of zeros. 
Another possibility is to use geometric methods to investigate the structures generated by 
the integral (in this way [T] allowed to find a new integrable system in [2]). 

Anyway the advantages of having an integral are so huge that the search for additional 
integrals of axially symmetric Ricci-flat metrics should be continued. But as our results 
indicate if the spacetime fll.ip is integrable, it should be rather a nontrivial phenomenon. 
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